The dynamics of quantum vortices in a toroidal trap 
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The dynamics of quantum vortices in a two-dimensional annular condensate are considered by 
numerically simulating the Gross-Pitaevskii equation. Families of solitary wave sequences are re- 
ported, both without and with a persistent flow, for various values of interaction strength. It is 
shown that in the toroidal geometry the dispersion curve of solutions is much richer than in the 
cases of a semi-infinite channel or uniform condensate studied previously. In particular, the toroidal 
condensate is found to have states of single vortices at the same position and circulation that move 
with different velocities. The stability of the solitary wave sequences for the annular condensate 
without a persistent flow are also investigated by numerically evolving the solutions in time. In 
addition, the interaction of vortex-vortex pairs and vortex-antivortex pairs is considered and it is 
demonstrated that the collisions are either elastic or inelastic depending on the magnitude of the an- 
gular velocity. The similarities and differences between numerically simulating the Gross-Pitaevskii 
equation and using a point vortex model for these collisions are elucidated. 

PACS numbers: ?? 



I. INTRODUCTION 

A Bose-Einstein condensate (BEG) offers an excellent 
quantum system to investigate, both experimentally and 
theoretically. Experimentally, virtually all of the physi- 
cal parameters can be controlled and thus the dynamics 
of quantum systems on a macroscopic scale are readily 
observed. It is a challenge, theoretically, to devise rele- 
vant and accurate models to accompany the increasing 
number of physical phenomena currently being reported. 
Quantum vortices are one such nonlinear physical phe- 
nomena that have been observed in BEGs [H |21 131 IH [3 E] , 
and that have been investigated theoretically [3 [HE]- In 
a BEG, the vortices have quantized circulation and their 
interactions on an inhomogeneous density background 
have been studied extensively pUl [TT] [H] [H} Hi] . 

Recent experimental interest [13 US] has shifted to fo- 
cus on BEGs in toroidal trapping potentials which are 
now more experimentally accessible. A BEG in a toroidal 
trap will form a ring shaped condensate. The first exper- 
imental observation of spontaneous vortex formation as 
the alkali gas is cooled through the transition tempera- 
ture to create a BEG has been reported in [16] for both a 
toroidal and a harmonic trap. In previous studies vortex 
formation through the transition temperature could be 
attributed to a number of different factors such as stir- 
ring of the condensate by a laser beam. The conditions 
under which a ring condensate exists and the manipula- 
tion of a condensate in a toroidal trap has recently been 
considered by [T7] and the generation of solitary waves 
have been studied in |18| . Experimental interest in the 
dynamics of superfluids in a ring is not new. Donnelly 
& Fetter [T3] noted that, in superfluid ^He, vortices ap- 
pear in the ring at a particular angular velocity when 
vortices can first compensate for the difference in irrota- 
tional velocity between the inner and outer boundaries 



of the ring. Experimentally, rotation of the condensate 
will create instabilities around the edge of the conden- 
sate. As the rotation velocity is increased the instabil- 
ities continue to grow and form vortices (see [20l [21]). 
A three-dimensional toroidal trapping potential, which 
creates a ring shaped condensate, might provide the nec- 
essary energy requirements to ensure stationary solitary 
wave solutions exist. 

Solitary wave sequences have previously been consid- 
ered in both two- and three-dimensional homogeneous 
condensates [22J [23]) a three-dimensional cigar shaped 
condensate [2U [25J [23] and also a two-dimensional chan- 
nel condensate [10]. The toroidal geometry considered 
here ensures that there are no phase differences in the 
condensate. As a result, the dispersion curve of the soli- 
tary wave sequences is far richer than has previously been 
seen. 

In a recent experiment, the observation of persistent 
flow in a toroidal trap and the existence of vortices, has 
been reported in [13]. Persistent flow in a BEG is the 
persistent circulation facilitated by the frictionless flow 
of a superfluid system. A useful way to view persistent 
flow is to consider a single vortex pinned to the centre of 
the ring shaped condensate. The energy at the centre is a 
(local) minimum and as such it will cost too much energy 
for the vortex to drift away from the centre. The vortex 
creates a constant circulation around itself: a persistent 
flow. A persistent flow in a BEG is also referred to as 
a supercurrent. The properties of a persistent flow in a 
toroidal trap have been theoretically considered in [27] • 
It is thought that studies on persistent flows in a BEG 
held under a toroidal trap could lead to crucial insights 
into the features of the critical velocity in superfluid "^He 
and to the fundamental relationship between superfluid- 
ity and Bose-Einstein condensation [15] . 

A two-dimensional condensate with external potential 
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trap of the form used by [TS] is considered 

Vext{x,y,z) = ^TO (w^x^ + cj^y^) + 

Voexp{~2ix' + y^)/wl), (1) 

for a constant potential Vq, waist wq (the axial distance 
from the laser beam's narrowest point), mass m and har- 
monic frequencies components {cOx , '-^y} will create the 
required ring shaped condensate. Reduction to a two- 
dimensional condensate is equivalent to a tightly trapped, 
or confined, condensate in the axial coordinate. Such a 
confinement is experimentally achievable and reduction 
to two-dimensional condensates has been theoretically 
considered in [28]. Previous theoretical investigations 
pS', "321 1311 132] have concentrated solely on solitary wave 
solutions in the absence of a persistent flow, in particular 
vortex solutions, in an idealised two-dimensional annular 
shaped condensate where the density is uniform within 
the annulus and zero outside of the annulus (a 'boxlike' 
toroidal trap) . Their analysis uses estabhshed point vor- 
tex model techniques from classical hydrodynamics and 
has shown that the motion of vortices is similar to the 
motion of classical vortices in an annular shaped domain 
[33]. The trapping potential ([T]) creates an inhomoge- 
ncous condensate that is markedly different to a 'boxlike' 
condensate. 

This paper will focus solely on the dynamics of solitary 
wave solutions in a two-dimensional condensate held un- 
der a toroidal trapping potential. Section II will intro- 
duce the required mathematical formulation for the an- 
nular condensate without a persistent flow. The families 
of solitary wave sequences that exist are then reported 
in Sect. III. In Sect. IV, the stability of the solitary wave 
sequences are numerically investigated and in Sect. V 
the evolution and collisions of vortex-vortex pairs and 
vortex-antivortex pairs are considered. Section VI adds 
a persistent flow to the annular condensate and reports 
the families of solitary wave sequences. Finally, the paper 
ends with a conclusion. Sect. VII. 



II. FORMULATION 

The dynamics of an annular shaped condensate are 
accurately described by the time-dependent dimensional 
Gross-Pit aevskii (GP) equation in terms of the macro- 
scopic wave function tp = ip{r, 6, t), 

= '^^'^ " {Ev - 6) - C/olVP) ^, (2) 

for two-dimensional external potential trap VextifyO)- 
The two-dimensional coupling constant is Uq, the mass 
of a boson is m and Ey is the chemical potential of the 
system. The external potential trap is taken to be solely 



in the radial direction such that 

V,xt{r,0) = V,xt{r) 

= Voexp{-2rywi) + ^mLj\^, (3) 

for a constant potential Vq, waist wq and frequency lo 
(lu = LOx ~ ^y, — 0). For small r the exponen- 
tial term will dominate over the harmonic term and, for 
a sufficiently tuned Vq, a circular hole in the centre of 
the condensate, an inner boundary, is created where no 
atoms will be present. Conversely, at larger values of r 
the harmonic term will dominate causing the condensate 
to form an outer boundary. Note that the toroidal trap 
can also be modelled by an external potential trap of the 
form 

Fe,t(r) = ^mc^2 (^A^ + ar^^ , (4) 

for (T = — 1 as has been considered by many authors 
(see for example [21]) for oscillator length h±^ and rel- 
ative strength of the quartic term A. Furthermore, a 
toroidal trap can be created by using Eq. |4jwith cr = -1-1 
and placing the condensate into rapid rotation. Above 
a certain critical angular velocity of rotation, a central 
hole is created in the condensate resulting in a toroidal 
geometry (see [SH [35l [36j [37] ) . However only ([3| will be 
used throughout this paper. 

Equation Q can be non-dimensionalised according to 

ip ^ V, ^^—^^^ r^aj_r, (5) 

where the transverse oscillator length is aj^ = 
(/i/v^ww)^/^ and the number density of the groimd 
state hsn — gfi^ / {2mUo) where the non-dimensional cou- 
pling potential g has been introduced. Thus the non- 
dimensional GP equation describing the dynamics in the 
annulus is 

2i^ = ~V^^^{^l~V{r)-g\^b\^)^, (6) 

for chemical potential /i = \f2Ey j huj , coupling constant 
g — \f2nUf^lfujja}y^ and where the external potential trap 
is 

l^(r) = Aexp(-?2^2)^ 1^2^ 

for potential A — \f2Vi^lhijj and inverse waist / = 
(v^/i/mwwg)^/^. Equation ([H]) is a four parameter sys- 
tem {g, A, /i} and is subject to the normalisation 

/ \i;\^rdrde ^271 \tl;\^rdr^l, (8) 

where V is the entire spatial domain. 
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A good approximation to the ground state for a par- 
ticular choice of parameters, {g, A, I, /z}, can be garnered 
from the Thomas-Fermi (TF) approximation, given by 



0, 



1/2 



^i>Vir), 
otherwise. 



(9) 



The four parameter set {g,A,l,fi} can be reduced to 
three by substitution of the TF approximation ([9| into 
the normalisation condition (|8]), i.e. 



27r 



/i — Aexp{—Pr^) — 



1„2 

2' 



'dr = 1, 



(10) 



such that the value of the chemical potential /i (say) is 
determined by the choice of {g, A,l}. An example of the 
ground state for g — 500, A — 100 and I — 0.9 is shown 
in Fig.[l] Note that since the Thomas- Fermi profile Q is 
only an approximation to the ground state, a numerical 
procedure, described in Sect. Ill, is required to find the 
exact ground state -00 for a particular choice of {g, A, I}. 
The corresponding chemical potential is found, by the 
same numerical procedure, to be /i = 14.72. 




FIG. 1: (Color online) The ground states of a condensate acting 
under a toroidal potential trap ^ for g = 500, A = 100 and I = 0.9. 
Frame (a) is a contour plot of the exact ground state ■00 while (b) 
is a slice along constant 9 of the exact ground state (solid line) 
and the Thomas- Fermi approximation 1p^F ^ (dashed line). 

The energy functional and angular momentum for the 
inhomogeneous condensate are given respectively by [10 , 



|V^|2 + (V{r) ~ + ^|^|Vdrd^,ll) 



{r-ro)^drde. (12) 



The energy of the solitary waves can be determined from 
E = Ef — Eg where Eg is the energy of the ground state 
in the absence of any solitary waves and is found numer- 
ically by evaluating Ef for ip = ipQ. Alternatively, an 
approximate analytically expression for Eg is found by 



evaluating Ef for ip = ■0tf- 

If solitary wave solutions are found to exist it is entirely 
plausible that a multitude of possible solutions could ex- 
ist for various 6. A restriction to the type of solutions 
sought is therefore taken so that only solutions with den- 
sity minima along 6 = ±7r/2 will be considered. Despite 
this restriction, the essential dynamics of solitary waves 
in an annulus will still be reflected. 



III. SOLITARY WAVES 

Solitary wave solutions of Eq. (|6| that preserve their 
form and move with a constant angular velocity at con- 
stant r subject to the external potential ([t]) are sought. A 
typical example of such a solitary wave solution is a single 
vortex, which will move in the inhomogeneous condensate 
because of the velocity field created by the surface of the 
condensate. Other solitary waves that are sought are a 
vortex pair, dark (gray and black) solitons and finite am- 
plitude sound waves - rarefaction waves. Thus, Eq. ^ IS 
recast in the frame rotating with the solitary wave using 
9' = e-m so that d/dt -Vld/de'. Hence 
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W^^+{fi~V{r)-g\tpf)ij, (13) 



where 6' is replaced by 9 for convenience and is solved 
subject to boundary conditions 



^ir,9,t)^Mr) 



{0,00}. 



(14) 



In the above statement for the boundary conditions, the 
ground state wavefunction at r = {0, 00} is zero, i.e. 
ip{0) = ip{oo) = 0. 

The solitary wave solutions are found numerically by 
a Newton-Raphson iteration technique. In order to nu- 
merically model the shape of the condensate, a cut in the 
infinite two-dimensional domain is taken along 9 — and 
the domain is unfurled so that the new two-dimensional 
rectangular domain in (r, 9) occupies the upper half 
plane. This semi-infinite numerical domain is mapped 
by the transformation r = tan~^(Dr) to a finite grid 
(0, 7r/2) X (0,27r) where is a constant chosen to lie 
in the range D ^ 0.4 — 0.8. The boundary conditions 
at = and 9 = 2tt are taken to be periodic. The 
resulting equations are expressed in second-order finite- 
difference form. Taking 201 x 200 grid points in the finite 
domain, the discretised nonlinear equations are solved by 
a Newton-Raphson iteration procedure using a banded 
matrix linear solver based on the biconjugate gradient 
stabilized method. The accuracy of the obtained solu- 
tions is verified by evaluating various integral identities. 
Indeed, an alternative expression for the energy can be 
found by taking 9 —* b9 for constant b in the expres- 



sions for energy and angular momentum (111 and (12 1 
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and considering the variational relationship 
d 



db 



[E - np) 



0. 



(15) 



b=l 



The result is an alternative expression for the energy 



E : 





dip 


Jv r 


do 



drd9, 



(16) 



which can be used as a check on the numerical results. 

An initial ansatz for the numerical procedure for wave 
function -ip can be obtained from 



(17) 



where ipQ is the exact ground state wave function (ob- 
tained numerically from the TF approximation ip^p) and 
ipsoi is motivated by the 'Tsuzuki' dark soliton solution 
in homogeneous [35] and inhomogeneous [TUl US] conden- 
sates 

ipsoi — ici + C2tanh(c3rcos6'), (18) 
for constants Ci, C2 and C3 to be chosen. Notice that 



( 18 1 creates a soliton at cos 6* = to ensure that solitary 
wave solutions exist only along 9 = ±tt/2. An alternative 
ansatz which looks directly for vortex solutions can be 
found from 



tp = IpQTpv, 

where the wave function for the vortex at (rp, tt) is 
r cos 6 + i{r sin 6 ~ tq) 



Ipv = 



+ {rsmO- ro)2 -|- 2/^ 



(19) 



(20) 



Additional vortices can easily be included into ( |19[ ) by 
writing ipy — HiLiV'ii where N is the number of vor- 
tices, as required. The position of the vortex is detected 
numerically by finding zeros of the real and imaginary 
parts of the wave function. 

As a vortex moves close to the boundary of the conden- 
sate, where ip —* 0, it becomes impossible to distinguish 
between a rarefaction wave and a vortex solution. In 
what follows the solitary wave solution is referred to as 
a "rarefaction wave" when the numerics cannot resolve 
whether or not the real (imaginary) part of ip changes its 
sign along the polar angle (see also [10] for discussion). 

An asymptotic expression for the group angular veloc- 
ity is obtained from the variation tp — s- ip + Sip in (11) 



and ( 12 1 and allows a comparison with the results of the 



numerics to be made. On use of the GP equation (13) 
one obtains 



dE 
dp ' 



(21) 



with the derivatives taken along the direction of propa- 
gation of the solitary waves. To obtain an approximation 
to the angular velocity as a function of ro, the position 



of the vortex, Eq. (21 1 is rewritten as 
dE _ dE/dro 



n = 



dp dp/dro 



(22) 



for the energy and angular momentum given as in (16 1 



and (12), respectively. 

Of the three parameters contained within the system 
{g, A,l} (the value of fi is determined from the other 
three), the interaction strength g can be considered to 
be the dominant factor. The interaction strength directly 
determines the number of healing lengths that span the 
condensate and thus in turn directly determines the num- 
ber of solitary wave solutions, along a constant 9, that 
can be present in the condensate. A range of values of 
g have been simulated (30 < g < 1000) of which two 
typical examples will be detailed below. In each example 
the energy-angular momentum dispersion curve is sym- 
metric about p = 0. The effect of this symmetry is to 
simply consider solitary wave solutions that move with 
the opposite circulation. In order to aid the analysis 
throughout this section different shades of colors will be 
used. Shades of black will represent a condensate with 
one solitary wave solution present, shades of dark gray 
(red) will represent a condensate with two solitary wave 
solutions present and shades of light gray (green) will 
represent a condensate with either three or four solitary 
wave solutions present. 



A. Low Interaction Strength 

A condensate with a low interaction strength will first 
be considered with the values of the parameters taken as 
{g,A,l} = {150,40,1} with the corresponding value of 
the chemical potential /i = 8.59. The complete family 
of solitary wave solutions in the annulus are then found 
and the energy-angular momentum dispersion curve is 
shown in Fig. |2] Two distinct solutions are found to ex- 
ist. The first is a single vortex (thick solid black line) 
of positive circulation (p > 0) which exists for angular 
velocity -0.29 <n < 0.40 and along 61 = 7r/2 [39]. The 
vortex will transcribe circles around the origin for con- 
stant angular velocity and for constant radius tq. At 

= 0.38, the single vortex loses its circulation and be- 
comes a rarefaction wave which is present for increasing 
angular velocities until the termination angular velocity 
is reached (for the range of parameters chosen here, the 
termination angular velocity is c = 0.40). No solutions 
are present above the termination angular velocity. 

A plot of the angular velocity of the vortex against dis- 
tance, ro, from the centre of the condensate is shown in 
Fig. [3] The shape of the curve is qualitatively similar (in 
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the bulk of the condensate) to that observed in a 'box- 
like' toroidal condensate (see [29] and [31]). However, at 
the boundaries of the condensate there is no qualitative 
agreement: the curve of Fig. [3] reaches finite values for 
the angular velocity at the inner and outer boundaries. 
Notice that the vortex has zero angular velocity when 
ro « 1.61. 
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FIG. 2: (Color online) The energy-angular momentum disper sion 
curve for the annular condensate governed by the GP equation | |13| l 
for parameter set {g,A,l} = {150,40, 1}. A single vortex solution 
for = 7r/2 is represented by the thick solid black line while the 
thick dashed black line represents a rarefaction wave for 6 = 7r/2. 
The thick dashed dark gray (red) line corresponds to a rarefaction 
wave for 9 = it7r/2. The dispersion curve is symmetric about p = 0. 



where now the angular velocity is increased, a second 
solution is found to exist as two rarefaction waves along 
9 = ±7r/2 respectively (thick dashed dark gray (red) line 
in Fig. [2]). The two rarefaction wave solution exists all 
the way to p = where a maximum on the dispersion 
curve is attained and the solution corresponds to a black 
soliton. 

A sample of the two distinct solutions on the dispersion 
curve is provided by the contour plots in Fig. [4] The 
single vortex of positive circulation at = 0.1 (frame 
(a)) and the two rarefaction waves at = 0.2 (frame (b)) 



are depicted. Note that by (21 1, the angular velocity is 



the gradient of the energy- angular momentum dispersion 
curve. 
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FIG. 4: (Color online) Two contour plots depicting two solutions 
of the GP equation ( |13[ l for parameter set {g,A,l} = {150,40, 1}. 
Frame (a) contains a single vortex at O = 0.1 (ro = 2.49) and frame 
(b) pair of rarefaction waves along 6 = ±7r/2 for Q = 0.2. 
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FIG. 3: The angular velocity Q plotted against the distance ro of 
the single vortex of positive circulation depicted in Fig.[2|for p > 0. 
The inner edge of the condensate is at r = i?i = 0.78 and the outer 
edge of the condensate is at r = R2 = 4.65. 

For < —0.29 the single vortex solution does not exist. 
Instead a cusp is formed similar to that observed in semi- 
infinite cigar-shaped condensates [TUlllS]. Over the cusp. 



B. Mid-Range Interaction Strength 

As the interaction strength is increased new features 
are seen to develop in the energy-angular momentum 
dispersion curve that are unexpected. An example of 
the features observed to occur can be seen by taking the 
parameter set {g,A,l} — {500,100,0.9} with the corre- 
sponding value of chemical potential fi = 14.72. The dis- 
persion curve is given in Fig. [5] and a selection of contour 
plots for the different solutions of the dispersion curve are 
given in Fig.'s [6] and [7j Note that the dynamics within 
this particular system are much more diverse than for the 
low interaction strength g — 150 detailed in the previous 
section. The reason for the increase in the number and 
type of solutions realised is the increase in the conden- 
sates' spatial extent. It is now possible for the conden- 
sate to contain at least two solitary wave solutions for 
any = ±tt/2. 

The dispersion curve contains seven distinct branches 
with each branch allocated a label (i)-(vii). At zero en- 
ergy and zero angular momentum there is no solitary 
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FIG. 5: (Color online) The energy-angular momentum dispersion 
curve for the annular condensate governed by the GP equation | |13^ 
for parameter set {g,A,Vj = {500, 100,0.9}. There are seven dis- 
tinct branches (i)-(vii) (the others can be obtained by symmetry 
about p = 0) . Branch (i) corresponds to a single vortex of positive 
circulation (thick solid black line), branch (ii) corresponds to two 
vortices on 6 = 7r/2 (thin solid dark gray (red) line), branch (iii) 
corresponds to a single vortex of negative circulation (thin solid 
black line), branch (iv) corresponds to two vortices at symmet- 
ric distances on 9 = ±-k/2 respectively (thin dashed-dot dark gray 
(red) line), branch (v) corresponds to the four vortex solution (thick 
solid light gray (green) line) , branch (vi) to a three vortex solution 
(thin dashed-dot light gray (green) line) and branch (vii) corre- 
sponds to two vortices at antisymmetric distances on 9 = ±-k/2 
respectively (thick solid dark gray (red) line). 



wave solution, however an increase in energy and (posi- 
tive) angular momentum will create a single solitary wave 
solution moving near the termination angular velocity at 
the far edge of the condensate along = 7r/2. The soli- 
tary wave solution is a single vortex of positive circula- 
tion. As the angular velocity is decreased, branch (i) in 
the dispersion curve is transcribed and the position of 
the vortex is seen to move towards the centre of the con- 
densate. At Q. ~ —0.11 a kink develops in the single vor- 
tex solution branch with energy and angular momentum 
given by {E,p) — (0.13,0.70). The same single vortex 
of positive circulation solution is found to be continuous 
over the kink. The solution is present until f2 = —0.12 
at which point the branch terminates with the vortex at 
ro = 2.35. 

A new solution is found to exist that joins with branch 
(i). The new branch, (ii), is characterised by the creation 
of a new vortex of negative circulation at the far edge of 
the condensate on 9 — 7r/2; see Fig. |6] frame (ii) for a 
contour plot of the new solution. Branch (ii) exists for 
—0.12 < < 0.13 and includes a solution where = 0. 
As the angular velocity increases both the vortex of pos- 
itive circulation (near the inner edge) and the vortex of 
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FIG. 6: (Color online) Contour plots of the three different config- 
urations (i)-(iii) given on the dispersion curve (Fig.[5]l for angular 
velocities fl = 0.1, S7 = and = 0.1 respectively. 

negative circulation (near the outer edge) move towards 
the origin. At 51 = 0.13, where the branch terminates, 
the positive vortex reaches the inner edge of the conden- 
sate and is lost while the negative vortex is at tq = 2.77. 

A new cusp is formed on the dispersion curve at 
{E,p) — (0.13, 0.34). As the angular velocity is decreased 
a new branch is traced out (branch (iii)) which is char- 
acterised by the existence of a single vortex of negative 
circulation. Further decreasing the angular velocity shifts 
the single vortex towards the outer edge of the conden- 
sate until a,t — —0.30 the vortex is lost and the branch 
terminates. 

A further solution can be found by increasing the an- 
gular velocity and creating a new branch, (iv). Branch 
(iv) is characterised by the appearance of two identical 
vortices at the far edge of the condensate on 9 = ±7r/2 
respectively; see Fig. [T] frame (iv) . The cusp formed by 
branches (iii) and (iv) is unique to all the other cusps 
in the dispersion curve. At the cusp the dynamics con- 
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FIG. 7: (Color online) Contour plots of the four different configu- 
rations (iv)-(vii) given on the dispersion curve (Fig. |5} for angular 
velocities fl = 0.2, Q = 0.1, Q = 0.1 and Q = 0.1 respectively. 
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FIG. 8: (Color online) Contour plots for f7 = for branches (v) 
and (vi) on the dispersion curve (Fig. [sjl showing the occurrence 
of a black soliton. 



tain no solitary wave solutions, similar to the termina- 
tion point of branch (i) at the termination speed (p — 0). 
In the experimental procedures it is possible that dis- 
turbances caused outside the far edge of the condensate 
could theoretically manifest into forming solutions exhib- 
ited by branches (i), (iii) or (iv). 

The angular velocity range of existence of branch (iv) 



is— 0.30 <f2< 0.13 during which the energy and angular 
momentum change considerably. As the angular velocity 
is increased the two vortices move towards the inner edge 
of the condensate each keeping the same fixed distance 
from the origin. The solution can be viewed as a sym- 
metric pair of vortices on 9 = ±tt/2 (there is symmetry 
about y = 0). The solution is present to negative angular 
momentum and terminates when {E,p) = (0.23,-0.38) 
with the vortices at tq = 2.54. 

At this point the dispersion curve forms a three-way 
cusp. Choosing to follow the thick solid light gray (green) 
curve in Fig. |5] (branch (v)) and thus decreasing the an- 
gular velocity, results in two new vortices being created 
at the edge of the condensate on 6 = ±7r/2, respectively, 
to produce a four-vortex configuration; see Fig. [7] frame 

(v) . The two new vortices are of opposite circulation 
compared to the two existing vortices. As the angular ve- 
locity is decreased towards zero the four vortices all move 
towards the inner edge of the condensate (the positions of 
the vortices on 9 = tt /2 always remain symmetric to those 
on 6* = — 7r/2). When the angular velocity becomes zero, 
two black solitons occur; see Fig. [s] frame (v). For neg- 
ative velocities the dynamics are symmetric to those for 
positive velocities such that when the branch terminates 
a,t n = -0.13 for {E,p) = (0.23,0.38), the configuration 
is identical to that at {E,p) = (0.23,-0.38). 

Starting at {E,p) — (0.23,0.38) and instead choosing 
to follow the thin dashed-dot light gray (green) curve in 
Fig. [5] and thus tracing out branch (vi), the penultimate 
configuration is realised. Here the dynamics consists of 
three vortices, with the new vortex appearing at the edge 
of the condensate for = 7r/2 (say). Thus the configu- 
ration consists of two vortices on 9 = tt/2 and a sin- 
gle vortex on 9 — — 7r/2; see the contour plot of Fig. [t] 
frame (vi). Increasing the angular velocity will cause all 
the vortices to gradually move towards the centre of the 
condensate until at zero angular velocity when the two 
vortices on 6* = 7r/2 decay to form a black soliton while 
the single vortex on 6* = —n/2 remains; see Fig. [S] frame 

(vi) . For positive velocities the single vortex continues 
to move towards the centre of the condensate, but the 
other two vortices now move towards the outer edge of 
the condensate. As the branch terminates at = 0.13, 
the furthest of these two vortices reaches the outer edge 
of the condensate and disappears. Meanwhile the two 
remaining vortices (on opposite sides of the condensate) 
are at different radii. 

The final realisable configuration is mapped out by 
branch (vii); see Fig. [t] frame (vii). The two vortices 
can be described as an anti-symmetric pair since they are 
at different radii (ri — 2.61, r-z = 3.38 for the vortices 
on 9 = +7r/2 and 9 = —7t/2 respectively for = 0.1) 
and their dynamics behave differently. As the angular 
velocity is decreased, the vortex at ri will move towards 
the inner edge of the condensate while the vortex at r2 
will move towards the outer edge of the condensate. At 
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zero angular velocity ri = r2 (note that there is no black 
soliton here). 

C. The Effect of A and I 

This section will consider in more detail the typical ef- 
fect that changing the values of the parameters A and I 
has on the dynamics. As explained before, the interac- 
tion strength g is the parameter that is most responsible 
for the width of the condensate and is thus the crucial 
parameter determining the type of dynamics that can 
be exhibited in the condensate. As a result it is thus 
expected that a change in A and I will not change the 
number of solitary waves that can span the condensate 
at any one time (since changing A and I is expected to 
only have a minor effect on the width of the condensate, 
and thus the number of healing lengths that span the 
condensate) . 

Thus consider the potential trap ([t]). Taking the 
derivative with respect to r and searching for the mini- 
mum one obtains 

dV o no 

— = -2yl/Vexp(-/V2 = 
or 

= +1-^ [ln(2A/2)]'/', (23) 

where the positive root is chosen and 2AP > 1, both 
to ensure r > 0. Thus the value of A and I effectively 
shift the position of the minimum of the potential and 
hence shift the position of the maximum of the density 
of the condensate. The effect will be to alter the range 
of angular velocities for which a particular solitary wave 
solution exists, however it is expected that A and I will 
have only a minor effect on the width of the condensate, 
which is dominated by the harmonic term in ([t]). 

A comparison of the energy-angular momentum dis- 
persion curve for two different parameter sets is shown 
in Fig. [9] for g — 500 and parameter sets: (a) {A,l} = 
{50,0.9} and (b) {A,l} = {50,1.5}. The dispersion 
curves can then be compared to Fig. |5] where {A, 1} — 
{100,0.9}. The value of g for all parameter sets is the 
same. Parameter sets (a) and (b) vary respectively one 
of the parameters {A, 1} of Fig. [5]while holding the other 
fixed. Thus a direct comparison into how A and I effect 
the dynamics of the condensate can be achieved. 

As can be seen in Fig. [9] the effect of changing A is 
almost negligible (compare Fig. [5] and Fig. |9ja)). Per- 
haps the only noticeable difference is a slight increase in 
the maximum energy, Emax, as A is decreased. The in- 
crease in Emax is to be expected since decreasing A will 
increase the width of the condensate slightly. The effect 
of changing I is again fairly minimal (compare frames (a) 
and (b) of Fig. [9]). However there are a couple of points 
that can be made. Firstly, increasing / has the effect of 
increasing Emax fairly considerably (from Emax = 0.28 
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FIG. 9: (Color online) A comparison of the energy-angular mo- 
mentum dispersion curve for two different parameter sets with 
g = 500. Frame (a) has {A,l} = {50,0.9} and frame (b) has 
{A, 1} = {50, 1.5}. The colors arc as in Fig. [5] 

on frame (a) to Emax — 0.35 on frame (b)). Secondly 
new features have developed in the dispersion curve for 
the three (thin dashed-dot light gray (green) curves) and 
four (thick solid light gray (green) curves) vortex solu- 
tions. Ignoring mirror symmetry about p ^ 0, the dis- 
persion curve in frame (a) contains a single three vortex 
configuration and a single four vortex configuration. In 
contrast the dispersion curve in frame (b) now contains 
two distinct three vortex configurations and four distinct 
four vortex configurations (see the blow-up of the region 
near Emax on Fig. [9](b)). Again the new configurations 
are to be expected; changing I will shift the position of 
the maximum density of the condensate without altering 
the width of the condensate much. 

The additional configurations that occur by changing I 
are important, but the essential dynamics of the problem 
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are unaltered; merely a few new configurations to already 
existing solutions are found (for instance the effect of in- 
creasing g might be to create new solutions where greater 
than four vortices in the condensate are present) . Indeed 
as (7 is raised to high interaction strengths many more 
distinct solutions will be present. 



D. Discussion 

The dynamics of a condensate held under a toroidal 
trapping potential have been described for two par- 
ticular cases: {g,A,l} — {150,40,1} and {g,A,l} — 
{500,100,0.9}. For {g,A,l} = {500,100,0.9} the solu- 
tions on branch (i) and branch (iii) have both been la- 
belled vortex solutions (see Fig. |5]) with the vortex on 
branch (i) possessing positive circulation and the vortex 
on branch (iii) possessing negative circulation. The con- 
tour plots of the two branches (see Fig. [6] (i) and (iii) re- 
spectively) indicate that the solutions are vortices. How- 



ever this raises an interesting anomaly. Figure 10 shows 
the angular velocity against distance plot for branches 
(i) and (iii). There is one point on the figure in which 
the curves intersect. At this intersection, the values of 
the angular velocity and the distances from the center 
of the condensate are equal for both branches, however 
the vortices possess different circulations. Such a sce- 
nario is mildly paradoxical and is certainly contrary to 
traditional vortex dynamical ideas. 




ro 



FIG. 10: The angular velocity Q plotted against the distance ro 
of the solutions of branches (i) (thick line) and (iii) (thin line) 
depicted in Fig. |5] for p > 0. The inner edge of the condensate 
is at r = Ri = 1.22 and the outer edge of the condensate is at 
r = R2 = 5.78. 

A plot of the real and imaginary components of the 
wave function {ip 



u + iv), Fig. 11 for the two branches 
shows that they do indeed contain vortex solutions (with 
the plots of the wave function across the direction of 
propagation confirming that the vortices have opposite 
circulation). However, closer inspection of the plots 
around r — 7 shows that the solution of branch (i) de- 



pletes to the ground state at a larger value of r than 
the solution of branch (iii). This indicates that branch 
(iii) possibly contains both a vortex solution and another 
embedded solution. Recall that branch (iii) was created 
when one of the vortices of branch (ii) reached the inner 
edge of the condensate and decayed into sound waves (in 
fact the vortex that decays into sound waves is the vor- 
tex of branch (i)). Therefore, the further solitary wave 
solution contained in branch (iii) is the remnants of the 
rarefaction waves created when this vortex decayed. 

Thus the solution of branch (i) can indeed be called 
a (true) vortex. However the solution of branch (iii) is 
actually a vortex with rarefaction waves embedded in the 
condensate. The statement resolves the apparent para- 
dox discussed earlier. To highlight their differences but 
at the same time to stress their similarities, the vortex 
of branch (i) will be denoted 'type I' and the vortex and 
rarefaction waves of branch (iii) will be denoted 'type IV. 
Type II solutions will often be referred to as vortex so- 
lutions in order to emphasize the similarities with type I 
solutions. 

For an annular shaped condensate the method of im- 
ages will produce two infinite sequences of image vortices. 
To see this consider a single vortex of positive circula- 
tion in the condensate at position r = tq with the inner 
edge of the condensate at i?i and the outer edge at R2 
so that Ri < tq < i?2- Firstly consider the image of 
the vortex with the inner edge of the condensate. By 
the Milne-Thomson circle theorem [ID] the image vor- 
tex, named In, will be of negative circulation and will 
be at r — i?J/ro. The image vortex In will itself have 
an image vortex, of positive circulation, associated with 
the outer edge of the condensate. The new image will be 
at r = R\rQ/R\ and is named 1 02- The process can be 
repeated with /02 now producing an image vortex with 
the inner edge of the condensate and so on. 

Secondly consider the image of the vortex with the 
outer edge of the condensate. The image vortex, named 
loi, will be of negative circulation and be at position 
r — R^/rQ. The image vortex Jqi will have an associated 
image vortex with the inner edge of the condensate to 
produce a vortex of positive circulation, named 7^2, at 
r = R\ro/R2- Again the process can be repeated. 

The upshot is that two infinite sequences of image vor- 
tices arc produced. For r < Ri, the first two terms of 
the sequence of image vortices are 



< 



(+) 



<^ <Ru 



(24) 



where a sign (±) has been added to indicate the circula- 
tion of the vortex. For R2 < r, the first two terms of the 
sequence of image vortices are 



> 



(+) 



Rj 



> 



(-) 



> i?2 



(25) 
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FIG. 11: Plots of the real (solid line) and imaginary (dashed line) 
components of the wave function (a) across the direction of prop- 
agation and (b) along the direction of propagation of the vortex 
for {g,A,l} = {500,100,0.9} and for the two single vortex cases 
type I and type II (branches (i) and (iii) in Fig. [s] respectively) for 
n = 0.1. 



It was seen in [IT] that the effect of a boundary in an 
inhomogeneous condensate was to alter the position of 
the image vortex by an amount equal to the depletion in 
density. The result was used in [10, for a channel con- 
densate to find an estimate for the velocity of the vortex 
as a function of distance from the centre of a channel. 
For the channel condensate there were again two infinite 
sequences of image vortices. However two distinct vor- 
tex solutions were not observed to occur. In the annular 
shaped condensate there is an asymmetry of the image 
vortices. The fact that the inner edge of the condensate is 
depleted by an exponential term and the outer edge of the 
condensate is depleted by a quadratic term distorts the 
positions of the image vortices and results in there being 
two distinct single vortex solutions. Of course there are 
other solutions that exhibit different configurations, for 
example {g,A,l} = {500,50,1.5} detailed in Sect. IIIC 
has four distinct four vortex solutions (see Fig. [9](b)). 

The geometry of the problem unfortunately does not 
lend itself to easy calculation of an analytical expres- 
sion for the angular velocity of a vortex as a function 
of distance r from the centre of the annulus. Various 
simplifications can be made in the hope of attaining an 
expression, among them being to consider only the two 
closest image vortices {In and loi)- However it has been 
out of reach to find an analytical expression for the an- 
gular velocity. 



IV. STABILITY 

The stability of the solitary wave sequences raises an 
intriguing question. Because quantized vorticity is an im- 
portant feature of condensate systems, the manipulation 
and observation of vortices is a crucial feature of any ex- 



periment involving vorticity. Here, a numerical approach 
is taken to investigate the stability of the solitary wave 
sequences in Sect. III. 

Families of solitary wave sequences were found in Sect. 
Ill for different values of the parameters {g,A,l}. The 
stability analysis in this section will concentrate on the 
parameter set {g,A,l} = {500,100,0.9} which was de- 
scribed in detail in Sect. IIIB. This parameter set is a 
typical example of the dynamics found in an annular 
condensate and the energy-angular momentum disper- 
sion curve is displayed in Fig. [5j 

To elucidate the stability of the seven distinct 
branches, a time evolution of the found solutions 
is performed. Specifically, the time-dependent non- 
dimensional Gross-Pitaevskii equation describing the dy- 
namics in the annular condensate ^ with toroidal trap 
([t]) is solved using a 4th order finite difference and 4th or- 
der Runge-Kutta scheme. The stationary solutions found 
in Sect. IIIB were found on an irregular polar grid. To 
translate the solutions to a regular cartesian grid, cubic 
splines are used. The solutions in each branch, for vari- 
ous angular velocities, are evolved in time until t ~ 150, 
a long enough time period to ensure a number of full 
revolutions of the solutions. It is expected that, over 
the course of a few revolutions, the unstable solutions 
will show one of two features: either they will decay into 
sound waves or, in the case of vortices will drift towards 
the inner or outer edge of the condensate. The stable 
solutions are expected to be unchanged after the time 
evolution. It is emphasised that the approach taken here 
does not prove the stability, or instability, of the respec- 
tive branches of the energy-angular momentum disper- 
sion curve, but this approach gives a strong indication to 
the respective stability of the different solutions. 

Of the seven branches of Fig. |5] it is observed that 
four branches are stable; branches (i) , (iii) , (iv) and (vii) , 
while three branches are unstable; branches (ii), (v) and 
(vi). The relative stability of the different branches is 
summarised in Fig. [12] The motion of the single vortex of 
type I and II (branches (i) and (iii) of Fig. [s] respectively) 
are both stable. A plot of the position of the vortex in 
each case is shown in Fig.[l3]for angular velocity = 0.1. 
One of the other stable branches, branch (iv), will be 
considered in more detail in Sect. VA. 

Turning now to the unstable solutions, branch (vi), 
where there are four vortices of alternate circulation 
along 9 = ±7r/2, is considered as an example of the 
dynamics observed. Figure [W] shows snapshots of the 
contour profile at different intervals of time for angular 
velocity = 0.1. As is seen, slight disturbances to the 
velocity profile cause the four vortex solution to decay. 
Initially, instabilities occur {t — 11), causing the four vor- 
tices to decay into two single vortices and sound waves 
(t = 12). After a short time interval, two vortices are 
observed to remain in the condensate which is now filled 
with sound waves (t = 26) . The decay of the vortices in 
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FIG. 12: The energy-angular momentum dispersion curve for the 
annular condensate for parameter set {g,A,l} = {500, 100,0.9} of 
Fig. [5] The relative stability of the different branches are indicated; 
stable (solid line), unstable (dashed line). 
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FIG. 14: (Color online) Contour profiles of the stationary solution 
for n = 0.1 of branch (vi) for four vortices placed along B = it7r/2 
at different intervals in time. The decay of the solution over time 
indicates instability. 

9 = 7r/2 and one vortex along 6 = — 7r/2, the instabil- 
ity develops along 9 = tt/2 while the single vortex along 
9 — — 7r/2 remains stable. However, over time the sound 
waves from the decay of the two vortices along = 7r/2 
propagate through the condensate and render the single 
vortex unstable. 



FIG. 13: The path (gray lines) of the single vortex solutions of 
type I and type II for n = 0.1 for integration forward in time. The 
vortices are at r = 3.11 and r = 2.87 respectively and are expected 
to be stable. The edges of the condensate are the black lines and 
are at r = i?i = 1.22 and r = R2 = 5.78. 



branches (ii) and (v) is similar. 

The relative stabilities of the seven branches is to be 
expected. Vortex rings/pairs in an infinite homogeneous 
condensate were found to be stable [23l [42] as well as 
vortex rings/pairs in a semi- infinite channel condensate 
[25j . However multiple vortex rings/pairs are unstable in 
these geometries. For the case considered in this paper, 
the unstable solutions are ones for which there exists two 
(i.e. multiple) vortices along the same polar angle. The 
two vortices will decay into a single vortex and sound 
waves (see Fig. 14 1. In contrast, solutions that contain 
only a single vortex along a polar angle are stable. In- 



deed, for branch (vi) where there are two vortices along 



V. EVOLUTION AND COLLISIONS OF 
SOLITARY WAVES 

The evolution of vortices in multiply connected do- 
mains has received much attention over the years. The 
majority of the focus has been on modelling classical fluid 
dynamical situations of interest, for example the motion 
of a vortex around islands |43| or through a gap in a wall 
[m US]. Traditional hydrodynamical methods for solv- 
ing such situations revolve around a point vortex model 
which enables a reduction in the dimension of the sys- 
tem whilst capturing all the essential dynamics. Point 
vortex models can also be successfully applied to annu- 
lar condensates; see [3^1 [SHI HU [12] • However the point 
vortex model relies on incompressible dynamics so the 
density in the bulk of the condensate is constant and is 
zero outside of the condensate. Thus, interactions be- 
tween the background condensate and the vortices are 
forbidden and so the transfer of energy between the vor- 
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tices and the background condensate does not occur. In- 
stead, vortices initially in the condensate will remain in 
the condensate. Decay of vortices is therefore not possi- 
ble using a 'boxlike' trap and a point vortex model. It 
has already been seen (Sect. IV) that interactions with 
the background condensate can have a significant effect 



posite sides of the condensate. As expected, the vortices 
traverse paths of constant radius and move with the same 
angular velocity. The paths of the vortices are shown in 



on the dynamics. Figure 14 shows that the interactions 
between the background condensate and the four vortices 
results in an instability developing and the decay of two 
of the vortices into sound waves. 

Experimentally, collisions of three-dimensional solitary 
waves have been observed in gB] and a recent theoret- 
ical study j37j considered the head on collisions of vor- 
tex rings and solitons in a cigar shaped condensate. It 
was observed in [IT] that, depending on the speed of ap- 
proach, the interactions of the solitary waves can be elas- 
tic or inelastic, with the solitary waves either repelling, 
passing through or annihilating each other. Collisions 
where the solitary waves pass through each other (elastic 
collision) is also forbidden by point vortex models. 

In a recent paper, Li et al. |14j considered vortex- 
antivortex pairs in a two-dimensional harmonically 
trapped condensate. In such a condensate it becomes 
difficult to simulate vortex-vortex or vortex-antivortex 
collisions. Instead a toroidal geometry explicitly allows 
this possibility. This section will consider the evolution 
of vortex-vortex pairs and vortex-antivortex pairs in the 
annular condensate by numerically simulating the GP 
equation ^ by the same method employed in Sect. IV. 
It is noted that, by solving Eq. ^ with toroidal trap 
([7]) , the bulk condensate density is not constant and thus 
the simulations in this section will provide a compari- 
son with the point vortex models employed in previous 
papers Hi [Ml 131 131] • 

The parameters of the condensate are again chosen to 
be {g,AJ} = {500,100,0.9}. The evolution and coUi- 
sional dynamics of two vortices are the main focus and 
thus the initial states for the condensate will contain two 
vortices, one along 6 — 7r/2 and one along 9 = —n/2. 
As detailed before, there exist two distinct single vor- 
tex solutions denoted type I and type II (branch (i) and 
branch (iii) on Fig. [sj respectively). The solutions on 
these branches are isolated (note that it is not the vor- 
tex solution that is isolated, but the complete solution in 
the absence of the ground state). The initial state can 
then be formed from -0 — ipoipv+tpv-, for ground state ipo 
and vortex states ip^^ (20) for vortices along 9 = ±7r/2 



respectively. Different combinations and different circu- 
lations of the single vortex solutions of type I and type 
II are considered. 



A. Vortex- Vortex Pair 

Consider first the evolution of two vortices of type I 
of the same circulation and placed at equal radii on op- 



Fig. 15 'a) for vortices placed at tq — 3.82 where the an- 
gular velocity is = 0.2. Also shown in Fig. [T5|[b) is the 
identical situation but now for vortices of type II placed 
at radii tq =4.38 where the angular velocity is = — 0.2. 
Note that the construction of these models forms conden- 
sates that are identical to branch (vi) of the dispersion 
curve in Fig. [5] which was shown to be numerically sta- 
ble. However, for a given radii of the two vortices, the 
angular velocity is different in all cases. 





FIG. 15: The path (gray Unes) of vortices placed initially on op- 
posite sides of the condensate with the same circulation. Frame (a) 
contains two vortices of type I of Fig. [s] at ro = 3.82 a,t 9 = it7r/2 
and frame (b) contains two vortices of type II at ro = 4.38. The 
edges of the condensate are the black lines and are at r = i?i = 1 . 22 
and r = R2 = 5.78. 



Next, the interactions of type I and type II vortices are 
considered. A vortex of type Iat6' = — 7r/2 and a vortex 
of type II at 9 — tt/2 are placed at the same radial posi- 
ton ro ~ 3.62. The vortices have the same circulation, 
however at rg = 3.62, the vortex of type I has positive 
angular velocity (fi = 0.18), whereas the vortex of type II 
has negative angular velocity (il = —0.1). Thus the vor- 
tices are expected to move towards one another but with 
different magnitudes of angular velocity. Indeed, as time 
evolves, the vortices begin to approach each other along 
the same radii with their initial angular velocities. How- 
ever as the distance between the vortices decreases, the 
vortices begin to lose angular velocity and move towards 
the centre of the condensate. They continue to move to- 
wards the centre of the condensate until they reach the 
inner boundary where they both decay into sound waves. 
A selection of snaps of the contour profile in Fig. [16] at 
different time intervals shows the process. 
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FIG. 16: (Color online) Contour profiles of the evolution of one 
vortex of type I placed at 9 = — 7r/2 of Fig. [sjand one vortex of 
type II at S = 7r/2. Both vortices are initially are at r = 3.62 and 
have the same circulation. 



B. Vortex- Antivortex Pair 

The collision of the two vortices of type I and type II, 
detailed in the previous section and in Fig.[T6]has opened 
up the question as to whether there is an analogy with 
the coUisional properties outlined by Komineas & Brand 
[47] , For solitary waves in a cigar shaped trap, the head- 
on collisions were elastic for low and high velocities of 
approach but were inelastic for intermediate values of ve- 
locity. To investigate the elasticity of the collisions in an 
annular condensate, a series of simulations for different 
angular velocities were considered for two vortices both 
taken to be of type I and placed on opposite sides of the 
condensate at the same radius but with opposite circu- 
lations. The three angular velocities chosen are = 0.1 
(with corresponding ro = 3.11), ft — 0.2 (rg = 3.82) and 
O = 0.3 (ro = 5.15). The three simulations are outlined 
in Fig. 's [17] [18] and [19] respectively. 

When the angular velocity is small (il — 0.1), initially 
the vortices approach one another with constant angular 
velocity. However as the distance between the vortices 
begins to reduce, their influence on each other begins to 
take effect and dominate the flow in the condensate. At 
such time, the angular velocity is reduced and the vor- 
tices move to a new radius where they begin to separate 
again at a new, constant, angular velocity. The process is 
repeated when the vortices next meet at the opposite side 



of the condensate, however now the vortices move back 
to the original radius and angular velocity. The process 
continues forever and is therefore elastic: the vortices 
upon collision merely transfer to a new energy and an- 
gular momentum configuration. In Fig. |17| a plot of the 
paths of the vortices is provided. 
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FIG. 17: The path (gray lines) of vortices of type I of Fig. [s] 
placed initially on opposite sides of the condensate with opposite 
circulation and with radius ro = 3.11. The edges of the condensate 
are the black lines and are at r = iJi = 1.22 and r = R2 = 5.78. 
The angular velocity is f! = 0.1. 



As the angular velocity is increased to f2 = 0.2, the 
vortices again upon collision slow down and move to a 
smaller radius. However, the vortices now possess too 
much energy to simply find a new radius to move along. 
Instead the vortices propel each other towards the in- 
ner edge of the condensate where they decay into sound 
waves. A similar process was observed to occur for the 
vortex-vortex pair considered in Fig. 



16 Contour plots 



in Fig. 18 at different time intervals for — 0.2 show the 
inelastic collision and decay process. 

Further increasing the angular velocity will cause the 
collisions to become elastic again. Taking = 0.3 as a 
typical example of the dynamics observed, the vortices 
now pass through one another and continue in their orig- 
inal direction of motion. However, during each collision 
a fraction of the energy of the vortices is lost to sound 
waves. Over time, therefore, the vortices slowly begin to 
lose energy. The energy lost during the collision to sound 
waves can be determined by 



AE = 



Ea — Eh 
Ea 



(26) 



where Ea is the energy before collision and £'{, is the 
energy after collision and where Eb is inferred from the 
angular velocity of the vortices after the collision. The 
value of the energy loss has been calculated after a num- 
ber of collisions for H. ~ 0.3 and is found to be always 



less than 1%. Figure 19 shows contour plots of the initial 
collision (between t = 5 and t — 6) and a plot of the state 
of the system after a long evolution (< = 150). 

The correspondence of the dynamics with the dynam- 
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FIG. 18: (Color online) Contour profiles of the evolution of vortices 
of type I of Fig. |5] Both vortices are initially are at ro = 3.82 and 
have opposite circulation. The angular velocity is f! = 0.2. 



FIG. 19: (Color online) Contour profiles of the evolution of vortices 
of type I of Fig.|5] Both vortices are initially arc at ro = 5.15 and 
have opposite circulation. The angular velocity is f2 = 0.3. The 
vortices pass through each other between t = 5 and t = 6. 



ics observed in |^ for the cigar shaped condensate is ev- 
ident. For the annular condensate, high and low angular 
velocities possess elastic collisions, whereas intermediate 
values of angular velocity possess inelastic collisions. 

Another striking feature of the interaction of the dy- 
namics in the annular condensate is the similarity of 
Fig.'s 



13 15 and 17 to the dynamics observed in 



by a point vortex model. Here, the figures show the mo- 
tion of the vortices as time is evolved. In inhomogeneous 
condensates where the density varies over all space, the 
motion of the vortices is because of the density gradi- 
ents' depletion at the boundary (see [TD]). The model 
detailed in Sect. IV and this section to find the evolution 
of vortices considers the depletion of the density gradi- 
ent. However, the point vortex model does not take into 
account the density gradient. The two models give the 
same trajectories but different velocities for each trajec- 
tory for both type I and II solutions. The point vortex 
model therefore provides a qualitative but not quantita- 
tive representation of the dynamics of the elastic colli- 
sions. Furthermore the point vortex model will not be 
able to capture any of the inelastic collisions. 



VI. PERSISTENT FLOW 

A persistent fiow (supercurrent) has been observed to 
occur in a toroidal trap by |TS], where the creation of 



vortices was also seen. The dynamics of solitary waves 
in a toroidal trap with a persistent flow present are read- 
ily obtained by a simple adaptation of the theory pre- 
sented in Sect. II. As noted before, a persistent flow can 
be characterised by the existence of a vortex at the ori- 
gin. Therefore take the non-dimensional GP equation ^ 
with the same trapping potential ([T]). Now consider the 
transformation of the wave function — tp exp{i0) which 
is equivalent to placing a singly quantized vortex at the 
origin. Then, the modifled GP equation now reads 



2i 



'dt 



2i dip 
^~d0 



1 



V'. 



Transfer to the rotating frame using 6' = 9 ^ 
d/dt 



(27) 
fii, so that 



-Vtd/dO' , and such that (27) becomes 



2i 



(hp_ 



n+\] =V',jj+[fi-V{r)-g\^\ 



(28) 

where the hat (") notation for ip and the dash (') notation 
for 6 has been dropped for convenience of notation. From 
now on ip will refer to the wave function in the condensate 



with a persistent flow. Equation (28) describes the dy- 



namics of solitary wave solutions in a toroidal trap with 
a persistent flow. Solitary wave solutions that move with 
a constant angular velocity f2 at constant r are sought by 
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the same numerical procedure introduced in Sect. III. 
The normahsation condition ([8| is unchanged; 



IV'pr drd6l = 27r / IV'prdr = 1, (29) 
V Jv 



where V is the entire spatial domain. The Thomas- Fermi 
approximation is now updated to be given by 



■0tf(?') 



^ — Acxp(— Pr^) — \ 



9 



1/2 



(30) 



while the modified energy functional is 



E, = i^lVVf +(v(r)-„+l)|Vf + 



which can be obtained from the transformation ip — > 
ijjexp(i9) in (111. The energy is then obtained from 
E = Ef — Eg, where Eg is the ground state energy and 



^9' 

is evaluated as En 



E 



Note that as expected 



dt Si;* ■ 



(32) 



The form of the angular momentum, however, is not 
so clear. Applying the transformation -ijj exp(i0) to 



( 12 1 will give 



Pi 



(^ - ^o) 



dip* 



(33) 

for ground state V'o- However this particular form of 
the angular momentum will be rejected and instead the 



following form (note its similarity to ( 12 1) of angular mo- 
mentum will be used 



P = 



(^ - V'o) 



d^* 
~d9 



(34) 



The reason for the rejection of ( 33 1 is to ensure that 



the angular momentum of the elementary excitations are 
found and not the angular momentum of the elementary 
excitations plus the angular momentum caused by the 
persistent flow. Additionally it allows the group rela- 
tionship = dE/dp to be satisfied for energy given by 
(31 1 and angular momentum given by (34 1. Finally, an 



alternative form for the energy can be obtained by tak- 
ing hO for constant b in the expressions for energy 



and angular momentum (31) and (34 1 and considering 



the variational relationship (21) to give 



E ■ 



1 

V r 



dip 
1)6 



2V 36 ^ dd 



drdB. (35) 



The dynamics of solitary wave solutions in a toroidal 



trap with a persistent flow have been simulated for a 
range of values of interaction strength 30 < g < 750 and 
parameters A and I. The results for two typical examples 
are detailed below for g = 30 and g = 150. As before, in 
each case the energy-angular momentum dispersion curve 
is symmetric about p = Q. The same colour scheme used 
in Sect. Ill will again be employed here. 



A. Interaction strength g = 30 

The first case study will focus on a parameter range 
with a low interaction strength: {5, A, 1} — {30, 15, 0.5} 
with the corresponding value of the chemical potential 
jjL — 8.24. The complete family of solitary wave solutions 
in the annulus are then found and the energy-angular 
momentum dispersion curve is shown in Fig. 



20 For the 



particular choice of parameters there is only one distinct 
solitary wave solution which exists in the range —0.08 < 
< 0.09 and is comprised solely of a single rarefaction 
wave along 6 = 1^/2 (see the thick black dashed line in 
Fig. 20). At no energy or angular momentum does the 



rarefaction wave develop vorticity. 



0.025 - 



-0.4 -0.3 -0.2 -0.1 



0.1 0.2 0.3 0.4 0.5 



FIG. 20: The energy-angular momentum dispersion curve for the 
annular condensate with persistent flow governed by the GP equa- 
tion ( |28| for parameter set {g,A,l} = {30,15,0.5}. Rarefaction 
waves are represented by the thick black dashed line. The disper- 
sion curve is symmetric about p = 0. 



B. Interaction strength g — 150 

The dispersion curve for {g, A, 1} — {30, 15, 0.5} was a 
typical example of the dynamics of the condensate in a 
persistent flow for low interaction strength. By increas- 
ing the interaction strength, a wider range of dynamical 
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solutions will be exhibited. The second choice of pa- 
rameter set is chosen to be {g,A,l} = {150,40,1} with 
chemical potential /i — 8.75. Note that this is the same 
parameter set as chosen in Sect. IIIA, however the chem- 
ical potentials are different. For the parameter set there 
exist three distinct solitary wave solutions which can be 
seen in the energy-angular momentum dispersion curve 
plotted in Fig. |2l] A plot of the angular velocity of the 
vortex solutions against distance from the origin is given 
in Fig. [22] 



E 




FIG. 21: (Color online) The energy-angular momentum dispersion 
curve for the annular condensate with persistent flow governed by 
the GP equation |28]l for parameter set {g,A,l} = {150,40,1}. 
There exist two distinct single vortex solutions: thick solid black 
and thin solid black lines with corresponding rarefaction waves 
(dashed thick black and thin black lines respectively). In addi- 
tion there is a two vortex solution (thin dashed-dot dark gray (red) 
line) with corresponding rarefaction wave (thin dashed dark gray 
(red) line). The dispersion curve is symmetric about p = 0. 



To describe the dynamics exhibited consider first the 
point {E,p) = (0,0) where no solutions exist. If the en- 
ergy and angular momentum are increased (to positive 
values) then a rarefaction wave develops along 6 — +7r/2 
and exists in the range of angular velocity 0.08 < < 
0.33. At = 0.08 the rarefaction wave develops vor- 
ticity such that the solution is now a single vortex of 
positive circulation. The single vortex solution exists un- 
til n = -0.31 at which point {E,p) = (0,0.68) and the 
solution terminates. The rarefaction wave and the single 
vortex solution just described form the first branch of the 
dispersion curve. 

A cusp develops at {E,p) — (0,0.68) and a new, sec- 
ond, branch is formed epitomised by the onset of two vor- 
tices of opposite circulation at the edge of the condensate, 
one along 9 = +7r/2 and the other along 6 — — 7r/2 creat- 
ing an anti-symmetric two vortex solution. The dynam- 




FIG. 22: (Color online) The angular velocity against distance of 
the three different vortex solutions of {g,A,l} = {150,40,1}. A 
slice along 6 = ±it/2 is taken. Colors refer to the different branches 
on the dispersion curve (Fig. |21|. 



branch that is formed exists for —0.31 < ft < 0.15 and 
is subdivided into two regions; the first region where the 
two vortex solution exists (—0.31 < < 0.06) and the 
second region where the two vortices lose circulation and 
become a pair of rarefaction waves (0.07 < fl < 0.15). 

The second branch terminates at = 0.15 at which 
point a cusp develops over which the dynamics are again 
continuous. Note that the continuity of the dynamics 
cannot be seen in Fig. [22] since at the cusp the solu- 
tions are rarefaction waves. Tracing out the new third 
branch in the dispersion curve gives a system that con- 
tains a single rarefaction wave that exists in the range 
0.04 < < 0.15. At = 0.04 the rarefaction wave de- 
velops vorticity and the solution is thus a single vortex 
which is present all the way until the branch terminates 
at {E,p) = (0,0). 

The dispersion curve for the parameter range chosen 
here therefore contains two distinct single solitary wave 



ics over the cusp are continuous (see Fig. 22 1 . The new 



solutions (thick black and thin black curves in Fig. 21 1 
as was previously seen in annular condensates with no 
persistent flow present (Sect. III). Note however that at 
{E,p) = (0, 0) where the branches of the two single soli- 
tary wave solutions meet the dynamics are not contin- 
uous, as can be seen in Fig. [21] where the first branch 
contains a rarefaction wave (thick black dashed line) and 
the third branch which contains a single vortex solution 
(thin black solid line). 

The two case studies {{g,A,l} — {30,15,0.5} and 
{150,40,1}) have concentrated on low-to-mid range in- 
teraction strengths. For higher g the energy-angular mo- 
mentum dispersion curves become too complicated to be 
described in a detailed yet concise fashion. However it 
can be appreciated that the net result of increasing g will 
be to produce new distinct solitary wave solutions and 
thus many branches in the dispersion curve. A further 
result of increasing g will be in the decrease in the range 
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of angular velocities for which the rarefaction waves exist. 



VII. CONCLUSION 

This paper has investigated the dynamics of a two- 
dimensional condensate held under a toroidal trapping 
potential. A toroidal trapping potential creates an an- 
nular condensate with richer and broader range of the 
solitary waves in comparison with the semi-infinite chan- 
nel geometries [TOl Hi] . 

The toroidal condensate is dependent on three inde- 
pendent parameters (that can be set experimentally) 
{g, A, I}, of which the interaction strength g is the dom- 
inant parameter. Changing the value of g determines 
the number of different types of solutions on the disper- 
sion curve. For low values of g there are two funda- 
mental branches of solutions (see Fig. [2] for {g, A, 1} — 
{150, 40, 1}). For a larger value of g, the number of fun- 
damental branches of solutions increases. The parame- 
ter set {g, A, 1} = {500, 100, 0.9} has been detailed ex- 
tensively in this paper. There are seven fundamental 
branches of solutions (see Fig. [5]) which contain up to 
four vortices in the condensate. An interesting scenario 
arises in which there exists two different branches each 
of which contains a single vortex with the vortices of 
opposite circulation. There is a solution where the two 
branches possess a vortex with the same angular velocity 
and are at the same distance from the origin. Since the 
vortices have opposite circulation, such a result seems 
rather surprising. One of the solutions was determined 
to contain not just the vortex, but in addition rarefac- 
tion waves embedded in the background of the conden- 
sate which resulted from the decay of a vortex from a 
different branch on the dispersion curve. 

The stability of the seven fundamental branches of 
the parameter set {500,100,0.9} was also investigated 
by evolving the solitary wave solutions forward in time. 
The solutions were evolved until three of four complete 
revolutions were completed, or the solutions decayed into 
sound waves. It was determined numerically that, if the 
solitary wave sequence contains more than one vortex 
for any given angle, then the branch is unstable. All the 
other branches were determined to be stable. 

By isolating the solitary waves from type I and type II 
solutions it is possible to investigate the possible range 
of collisional effects and compare them with a point vor- 
tex model. The solutions in this paper are all obtained 
by numerically simulating the Gross-Pitaevskii equation 
for the inhomogeneous condensate, whereas point vortex 
models consider a 'box-like' condensate where the density 
is uniform in the annular region and is zero outside this 
region. It is seen in this paper that collision of the vor- 
tices result in two distinct behaviours: either the vortices 
collide elastically or collide inelastically. Elastic collisions 
occur for low angular velocities where the vortices repel 



each other or for high angular velocities where the vor- 
tices pass through one another. For intermediate angular 
velocities, the vortices collide inelastically and decay into 
sound waves. Point vortex models are only able to cap- 
ture the collisions for low angular velocities. 

Finally, a separate condensate is considered where now 
a persistent flow is included into the condensate. It is in- 
teresting to note the similarities and differences between 
the annular condensate without a persistent flow (Sect. 
Ill) and one with a persistent flow (Sect. VI). 
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